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CU2Cdi_xZnxSnSe4  solid  solutions  were  synthesized,  and  their  phase  constitutions  and  thermoelectric 
properties  were  investigated.  The  solid  solutions  crystallized  in  the  stannite-type  structure  for  Zn  con¬ 
tents  x  up  to  0.65  and  in  the  kesterite-type  structure  for  0.7  <  x  <  1.0.  The  lattice  parameter  a  and  cell 
volume  V  of  the  compounds  decreased  linearly  with  increasing  x  for  both  the  stannite-type  (0  <  x  <  0.65) 
and  the  kesterite-type  (0.7  <  x  <  1)  structures.  The  lattice  parameter  c  decreased  with  increasing  x  for 
the  compounds  with  the  kesterite-type  structure  but  increased  for  the  compounds  with  the  stannite- 
type  structure.  The  c/a  ratio  increased  with  increasing  Zn  content,  which  indicated  an  weakening  of 
the  lattice  distortion.  The  Seebeck  coefficient  tended  to  decrease  with  increasing  Zn  content,  whereas  the 
electrical  conductivity  and  thermal  conductivity  increased.  The  figure  of  merit  ZT  increased  with 
increasing  x  over  the  composition  range  of  0  <  x  <  0.60  and  then  fluctuated  with  a  further  increase  in  x.  A 
maximum  ZT  of  0.23  was  achieved  for  Cu2Cdo.4Zno.6SnSe4  at  720  K. 
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1.  Introduction 

The  thermoelectric  effect  is  the  direct  conversion  of  a  temper¬ 
ature  difference  into  electrical  voltage  and  vice-versa  on  the  basis 
of  the  Seebeck  effect  or  Peltier  effect,  respectively.  Thermoelectric 
materials  can  be  used  as  power  generators,  coolers,  and  thermal 
sensors  [1,2],  The  efficiency  of  a  thermoelectric  device  depends  on 
the  dimensionless  figure  of  merit  ZT:  ZT  =  S2(jT/k  [3,4],  where  S  is 
the  Seebeck  coefficient,  k  is  the  thermal  conductivity,  a  is  the 
electrical  conductivity,  and  T  is  the  absolute  temperature.  Thus,  the 
thermoelectric  performance  of  a  material  can  be  improved  by 
enhancing  its  Seebeck  coefficient  and  electrical  conductivity  and  by 
reducing  its  thermal  conductivity  [5], 

High-performance  thermoelectric  materials  that  operate  in  in¬ 
termediate-temperature  ranges  often  contain  lead  (Pb),  especially 
p-type  thermoelectric  materials  such  as  PbTe— SrTe  [6]  (ZT  =  2.2  at 
915  K)  and  AgPbmSbTe2+m(LAST-rn)  (ZT  =  2.2  at  800  K)  [7],  It  is 
well-known  that  the  presence  of  lead  in  the  environment  causes 
toxic  effects  in  plants  and  animals.  Therefore,  exploring  novel 
thermoelectric  materials  with  high  ZTs  that  are  lead  free  and  that 
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are  composed  of  elements  that  are  abundant  in  the  Earth's  crust 
remains  an  interesting  topic  for  researchers. 

Cu-based  intermetallics  or  minerals  have  been  found  to  possess 
promising  thermoelectric  properties  [8-10].  The  quaternary  com¬ 
pounds  Cu2ZnSnSe4  and  Cu2CdSnSe4  have  aroused  considerable 
research  interest  from  scientists  in  recent  years  [11].  These  mate¬ 
rials  have  notable  thermoelectric  properties,  such  as  stabilized  p- 
type  conductivity,  a  relatively  high  Seebeck  coefficient,  and  a 
relatively  low  lattice  thermal  conductivity.  The  compounds 
Cu2ZnSnSe4  and  Cu2CdSnSe4  are  potential  intermediate- 
temperature  thermoelectric  materials.  The  ZT  values  reach  0.9  at 
860  I<  for  Cu2.iZno.gSnSe4  [12],  0.95  at  850  K  for  QQZnSno.9Ino.1Se4 
[13],  0.3  at  720  I<  for  Cu2.o5ZnSno.95Se4  [14],  and  0.65  at  723  K  for 
colloidal-synthesized  Cu2CdSnSe4  [15], 

Both  the  Cu2ZnSnSe4  and  Cu2CdSnSe4  quaternary  compounds 
possess  body-centered  tetragonal  symmetry.  The  Cu,  Cd,  Zn,  and  Sn 
atoms  are  all  coordinated  with  Se  atoms  in  a  tetrahedral  manner.  The 
oxidation  states  of  Cu,  Cd,  Zn,  Sn,  and  Se  in  the  structure 
are  +1,  +2,  +2,  +4,  and  -2,  respectively  [16],  The  compound 
Cu2ZnSnSe4  crystallizes  in  the  kesterite-type  crystal  structure  with 
the  space  group  14  as  confirmed  by  powder  X-ray  diffraction  (XRD) 
and  neutron  scattering  [17]  measurements,  whereas  Cu2CdSnSe4 
crystallizes  in  the  stannite-type  structure  belonging  to  space  group 
/42m  [18],  The  differences  in  atomic  sites  between  these  two  com¬ 
pounds  can  be  described  as  follows.  When  all  of  the  Sn  atoms 
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exchange  sites  with  half  of  the  Cu  atoms  in  the  Cu2CdSnSe4  /42m  cell, 
the  structure  transforms  from  the  stannite-type  to  the  kesterite-type. 
Cu2CdSnSe4  can  be  divided  into  two  structural/functional  units:  a 
Cu2Se4  tetrahedral  array  that  acts  as  an  electrically  conducting  unit 
and  a  CdSnSe4  tetrahedral  array  that  acts  as  an  insulating  unit, 
satisfying  the  requirements  for  being  classified  as  a  Slack's  phonon 
glass  electron  crystal  (PGEC)  [19],  PGEC  materials  exhibit  high  elec¬ 
trical  conductivity,  similar  to  a  crystal,  but  low  thermal  conductivity, 
similar  to  a  glass.  The  electrically  insulating  CdSnSe4  unit  responsible 
for  the  low  thermal  conductivity  is  inserted  between  two  electrically 
conducting  Cu2Se4  units,  which  account  for  the  electrical  conduc¬ 
tivity  [12],  Recently,  considerable  interest  has  arisen  in  the  study  of 
the  thermoelectric  properties  of  wide-band-gap  quaternary  com¬ 
pounds  because  conventional  thermoelectric  materials  exhibit  a 
narrow  band  gap,  where  a  bipolar  effect  may  reduce  the  thermo¬ 
electric  efficiency  [14],  The  reported  band  gap  Eg  is  1.44  eV  for 
Cu2ZnSnSe4  and  0.96  eV  for  Cu2CdSnSe4  [20], 

With  respect  to  the  thermoelectric  properties  of  the  compounds 
Cu2ZnSnSe4  and  Cu2CdSnSe4,  Cu2ZnSnSe4  exhibits  lower  resistivity, 
whereas  Cu2CdSnSe4  manifests  lower  thermal  conductivity.  The 
structural  transition  and  thermoelectric  properties  of  the  alloys  in 
the  Cu2ZnSnSe4-Cu2CdSnSe4  pseudo-binary  system  have  not  yet 
been  reported.  In  terms  of  element  doping  or  composite  materials, 
Cu2CdxZni_xSnSe4  solid  solutions  or  alloys  in  the  two-phase  region 
may  exhibit  good  thermoelectric  properties  combined  with  the 
advantages  of  both  the  end  compounds.  In  this  work,  we  prepared  a 
series  of  Cu2CdxZn1_xSnSe4  samples  with  x  =  0  to  1  using  melt¬ 
quenching  and  spark  plasma  sintering  (SPS).  The  crystal  structure 
and  thermoelectric  performance  of  the  samples  were  investigated 
in  detail. 

2.  Experimental  details 

High  purity  Cu  (99.99%),  Zn  (99.99%),  Cd  (99.99%),  Sn  (99.99%), 
and  Se  (99.99%)  were  used  as  the  starting  materials:  they  were 
weighed  in  the  desired  ratios  and  sealed  in  quartz  tubes  under 
vacuum  (<10-2  Pa).  The  tubes  were  placed  in  a  furnace  and  slowly 
heated  to  1273  K  at  a  heating  rate  of  1  K  min-1.  The  alloys  were  held 
at  1273  K  for  24  h  were  subsequently  cooled  slowly  to  773  K  and 
maintained  at  this  temperature  for  120  h  before  being  quenched  in 
liquid  nitrogen.  The  quenched  alloys  were  powdered  in  an  agate 
mortar  and  then  ball-milled  in  a  planetary  ball  mill  (QM-4F, 
Nanjing  University,  China)  at  200  rpm  for  12  h  using  a  hard  stainless 
steel  vial  and  zirconia  balls.  The  weight  ratio  of  balls  to  powder  was 
maintained  at  approximately  6:1,  and  the  mill  vial  was  evacuated 
and  then  filled  with  a  pure  H2  atmosphere  to  prevent  the  powder 
from  oxidizing  during  the  milling  process.  The  powders  were 
pressed  into  pellet  form  and  annealed  at  773  K  for  200  h  in  sealed, 
evacuated  quartz  tubes.  The  pellets  were  then  ground  into  a  fine 
powder  and  consolidated  by  SPS  at  923  K  for  5  min  under  a  vacuum 
of  10-2  Pa  or  better  and  a  pressure  of  48  MPa  in  a  graphite  die  with  a 
peak  current  impulse  value  of  703  A.  The  density  of  the  bulk  sample 
obtained  after  SPS  was  greater  than  97%  of  the  theoretical  value.  Bar 
specimens  with  typical  dimensions  of  13.0  mm  x  4.0  mm  x  4.0  mm 
were  prepared  for  electrical  property  measurements,  and  disk 
specimens  with  a  diameter  of  10.0  mm  and  a  thickness  of  2.0  mm 
were  prepared  for  thermal  conductivity  measurements. 

The  phases  in  the  samples  were  analyzed  by  X-ray  diffraction 
(XRD)  on  a  Bruker  D8  Advance  SS/18  kW  diffractometer  equipped 
with  a  Cu  Ret  radiation  source;  the  samples  were  scanned  in  the  28 
range  between  10  and  120°.  Rietveld  refinement  of  the  XRD  pat¬ 
terns  was  performed  using  the  software  Topas  3.1  [21  ].  The  Seebeck 
coefficient  S  and  the  electric  conductivity  a  were  measured  on  a 
ZEM-2  (Ulvac-Riko,  Japan)  under  a  helium  atmosphere;  the  mea¬ 
surement  errors  were  approximately  7%  and  10%,  respectively.  The 


thermal  diffusivity  (a)  and  the  specific  heat  capacity  Cp  were 
measured  simultaneously  via  a  laser  flash  method  on  a  TC-9000h 
thermal  constants  measuring  system  under  a  vacuum  of 
1.0  x  10-3  Pa.  The  bulk  density  pm  of  the  sample  was  calculated 
from  the  sample's  geometry  and  mass.  The  thermal  conductivity  k 
of  the  samples  was  calculated  according  to  the  equation  k  —  aCppm; 
the  estimated  error  of  the  measurements  was  approximately  13%. 
DSC  experiments  were  performed  on  a  SETARAM  16-18  Evolution 
using  a  constant  heat  rate  of  5  K  min-1. 

3.  Results  and  discussion 

3.1.  Crystal  structure 

The  XRD  patterns  of  the  SPS-sintered  Cu2Cdi_xZnxSnSe4 
(x  =  0—1.0)  samples  are  presented  in  Fig.  1.  These  patterns  can  be 
identified  as  resulting  from  the  compound  with  the  space  group  /- 
42m  or  1-4  in  the  tetragonal  system,  without  any  contribution  of  a 
second  phase  (Fig.  1(a)).  The  intensities  of  the  strongest  peak  (112) 
of  the  patterns  were  normalized  to  evaluate  the  intensity  of  the 
other  peaks  as  a  function  of  the  zinc  content  x.  Fig.  1(a)  reveals  that 
most  of  the  peaks  shifted  to  higher  angles  with  increasing  zinc 
content  x.  The  (220)  and  (024)  peaks  at  a  28  angle  of  approximately 
45.0°  and  the  (116)  and  (132)  peaks  at  a  28  angle  of  approximately 
52.5°  split  in  the  pattern  of  the  sample  with  x  =  0  but  overlapped 
gradually  as  x  was  increased  to  1.0.  With  increasing  zinc  content  x, 
the  intensities  of  the  (002),  (110),  (022),  and  (114)  peaks  decreased; 
however,  those  of  the  (011),  (013),  (121),  and  (123)  peaks  increased 
substantially  (Fig.  1(b)). 

We  cannot  distinguish  between  two  symmetries  of  1-4  and  I- 
42m  on  the  basis  of  the  XRD  results  because  the  extinction  con¬ 
ditions  for  these  two  space  groups  are  identical.  However,  we  can 
identify  them  by  Rietveld  refinement  on  the  basis  of  the  reliability 
factors  Rwp  and  Rfjragg,  which  are  calculated  from  the  intensity  using 
the  symmetry  and  atomic  occupations.  Rwp  reveals  the  goodness  of 
the  whole  profile  fit,  and  RBragg  provides  a  valuable  indication  of  the 
suitability  of  the  structural  model.  In  the  refinement  process,  Zn 
and  Cd  atoms  are  fixed  to  occupy  the  same  site;  however,  their 
occupancies  are  free  for  refining.  Both  the  stannite-  and  the 
kesterite-type  structure  models  were  used  to  refine  each  XRD 
pattern.  As  illustrated  in  Fig.  1(a),  the  patterns  were  divided  into 
two  groups:  group  A  for  the  samples  with  x  =  0,  0.1,  0.2,  0.3,  0.5, 
0.55,  0.6,  or  0.65  and  group  B  for  the  samples  with  x  =  0.7,  0.8,  0.9, 
or  1.0.  When  the  stannite-type  structure  model  for  the  samples  in 
group  A  or  the  kesterite-type  structure  model  for  the  samples  in 
group  B  was  used,  the  Rwp  was  approximately  3%,  the  RBragg  was 
0.35-1%  (Table  1),  and  the  occupancies  of  the  Zn  and  Cd  atoms 
were  close  to  the  nominal  composition.  In  contrast,  when  the 
kesterite-type  structure  model  was  used  for  the  samples  in  group  A 
or  the  stannite-type  structure  model  for  the  samples  in  group  B,  the 
Rwp  increased  to  more  than  5%,  the  RBragg  increased  to  more  than 
30%,  and  the  occupancies  of  the  Zn  and  Cd  atoms  exhibited  obvious 
deviations  from  the  nominal  composition.  Because  of  the  strong 
diffraction  intensity  collected  by  the  LynxEye  array  detector,  the 
goodness  of  fit  factor  (RWp/Rexp)  was  approximately  2—3,  which  is 
slightly  larger  than  2.  However,  the  low  values  of  Rwp  and  RBragg  and 
the  consistency  of  the  occupancies  of  the  Zn  and  Cd  atoms 
confirmed  the  reliability  of  the  results  obtained  by  Rietveld 
refinement.  The  results  confirmed  that  the  samples  in  group  A 
belong  to  the  stannite-type  structure  with  space  group  I- 42m, 
whereas  the  samples  in  group  B  are  of  the  kesterite-type  structure 
with  space  group  1-4  (Fig.  1(a)). 

The  FWHM  of  the  (112)  peak  with  the  maximum  intensity  ob¬ 
tained  by  fitting  of  the  peak  profile  shows  that  there  is  no  obvious 
broadening  of  the  (112)  peak  of  the  samples  with  different  x  values 
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Fig.  1.  (a)  Powder  XRD  patterns  of  G^Cdt-xZnxSnSe,!;  (b)  partial  details  of  the  XRD  patterns  (the  stacking  order  is  the  same  as  that  in  (a));  (c)  FWHM  of  the  (112)  peaks;  (d)  XRD 
Rietveld  refinement  results  for  the  x  =  0.6  alloy. 


Table  1 

Lattice  parameters  (a,  c,  and  V ),  c/a  ratio,  mean  atomic  weight  (M),  calculated  density  ( pc ,  g/cm3),  melting  point  (Tm),  value  of  (T^/2p2/3 ) /M7/6  (jt0,  R3,2g2/3  cm*2),  and  reliability 
factors  (Rexp,  Rwp,  and  RBragg,  %)  for  the  Rietveld  refinements  of  CujCd,  xZnxSnSe4. 


x 

o(A) 

c(A) 

V(A3) 

c/a 

M 

Pc 

Rexp 

R»p 

RBragg 

TJK) 

K0 

0 

5.829(1) 

11.403(1) 

387.5(1) 

1.956 

84.25 

5.776 

0.94 

2.40 

1.01 

1040 

612.04 

0.1 

5.813(1) 

11.409(1) 

385.5(1) 

1.962 

83.67 

5.764 

0.85 

2.27 

0.61 

1044 

619.76 

0.2 

5.797(1) 

11.415(1) 

383.6(1) 

1.969 

83.08 

5.753 

0.83 

2.44 

0.60 

1047 

626.78 

0.3 

5.782(1) 

11.418(1) 

381.7(1) 

1.975 

82.49 

5.741 

0.83 

2.36 

0.35 

1048 

632.02 

0.5 

5.756(1) 

11.421(1) 

378.5(1) 

1.984 

81.31 

5.718 

0.84 

2.38 

0.41 

1053 

645.57 

0.55 

5.745(1) 

11.422(1) 

377.0(1) 

1.988 

81.02 

5.616 

0.82 

2.93 

0.58 

1050 

637.83 

0.6 

5.740(1) 

11.419(1) 

376.2(1) 

1.989 

80.73 

5.604 

0.81 

2.78 

0.75 

1052 

641.46 

0.65 

5.732(1) 

11.414(1) 

375.1(1) 

1.991 

80.43 

5.587 

0.80 

2.82 

0.68 

1046 

637.40 

0.7 

5.727(1) 

11.410(1) 

374.2(1) 

1.992 

80.14 

5.686 

0.82 

2.89 

1.00 

1045 

646.74 

0.8 

5.714(1) 

11.392(1) 

372.0(1) 

1.993 

79.55 

5.676 

0.84 

2.29 

0.53 

1049 

655.30 

0.9 

5.703(1) 

11.370(1) 

369.8(1) 

1.994 

78.96 

5.666 

0.83 

233 

0.55 

1057 

667.78 

1 

5.693(1) 

11.345(1) 

367.7(1) 

1.993 

78.38 

5.657 

0.82 

2.99 

0.64 

1061 

676.74 
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(Fig.  1(c)).  This  result  indicates  that  all  of  the  samples  are  a  single 
stannite-  or  kesterite-type  phase.  The  two-phase  region  between 
the  stannite-  and  the  kesterite-type  solid  solutions  appears  narrow, 
between  x  =  0.65  and  x  =  0.70.  Notably,  however,  slight  atomic 
disorder  may  exist  in  the  samples  with  x  =  0.55,  0.60,  0.65,  or  0.70 
because  of  the  relatively  higher  Rwp  or  f?Bragg-  When  the  site  dis¬ 
order  was  refined  for  the  Cu  and  Sn  atoms,  approximately  12%  of 
the  Sn  atoms  exchanged  sites  with  Cu  atoms  in  the  sample  with 
x  =  0.7,  and  the  Reragg  factor  decreased  from  1%  to  0.76%.  Because 
the  atomic  scattering  factors  of  Cu  and  Zn  are  similar,  as  are  those  of 
Sn  and  Cd,  the  degree  of  disorder  cannot  be  determined  from  the 
XRD  Rietveld  refinement.  As  previously  described,  the  main 
structural  difference  between  the  stannite-type  Cu2CdSnSe4  and 
the  kesterite-type  Cu2ZnSnSe4  is  the  exchange  of  Sn  atoms  with 
half  of  the  Cu  atoms.  Therefore,  the  structural  transformation  be¬ 
tween  the  stannite-type  and  kesterite-type  may  be  accomplished 
using  atomic  disorder  in  Cu2Cdi_xZnxSnSe4  alloys. 

The  lattice  parameters  and  reliability  factors  for  the  compounds 
obtained  by  Rietveld  refinement  are  listed  in  Table  1,  and  the 
refined  atomic  parameters  for  Cu2CdSnSe4  and  Cu2ZnSnSe4  are 
listed  in  Table  2.  A  representative  XRD  Rietveld  refined  pattern  for 
the  sample  with  x  =  0.6  is  presented  in  Fig.  1(d).  The  lattice  pa¬ 
rameters  are  plotted  as  a  function  of  the  zinc  content  x  in  Fig.  2.  It  is 
observed  that  the  lattice  parameter  a  and  cell  volume  (V)  decrease 
linearly  with  increasing  Zn  content  x  for  the  compounds  in  both 
groups,  but  with  different  slopes.  This  result  illustrates  the  substi¬ 
tution  of  larger  Cd2+  ions  (0.78  A)  by  smaller  Zn2+  ions  (0.60  A)  in 
the  structure  [22],  The  lattice  parameter  c  increases  linearly  for  the 
stannite-type  compounds  but  decreases  for  the  kesterite-type 
compounds  with  increasing  x.  This  finding  may  account  for  the 
bonding  anisotropy  in  the  ab  plane  and  along  the  c  direction.  The 
linear  relationships  between  the  lattice  parameters  and  zinc  con¬ 
tent  x  can  be  expressed  as  follows:  as(A)  =  5.828(1)  -  0.149(3 )x; 
cs(A)  =  11.408(3)  +  0.0177(9)x  and  1/S(A3)  =  387.549(4)  -  19.098(1) 
x  for  stannite-type  compounds,  and  a«( A)  =  5.805(4)  -  0.113(4)x; 
CkH A)  =  11.564(8)  -  0.217(9)x  and  Vk( A3)  =  389.4(2)  -  21.7(3)x  for 
kesterite-type  compounds.  The  lattice  parameter  a  also  exhibits  a 
linear  relationship  with  the  mean  atomic  weight  (M)  of  the  com¬ 
pound  as  as(A)  =  0.025(1  )M  +  3.73(6);  a*( A)  =  0.019(1  )M  +  4.18(6), 
which  is  similar  to  the  formula  reported  by  Matsushita  et  al.  [20], 
where  the  subscript  S  refers  to  the  stannite-type  structure  and  K 
refers  to  the  kesterite-type  structure.  If  we  view  the  structure  of 
Cu2CdSnSe4  as  being  derived  from  the  cubic  zinc  blende  structure, 
the  lattice  distortion  can  be  estimated  from  the  c/a  ratio.  When  the 
zinc  content  x  increases  from  0  to  1.0,  the  c/a  ratio  gradually  ap¬ 
proaches  2  (Table  1 ),  which  means  that  a  larger  incorporation  of  Zn 
atoms  leads  to  less  lattice  distortion. 

The  atomic  distances  between  Se  and  M  (M  =  Cu,  Zn/Cd,  or  Sn) 
are  plotted  as  a  function  of  the  zinc  content  x  in  Fig.  3.  All  the  Se— M 
bond  lengths  decrease  with  increasing  x  in  both  the  stannite-  and 
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Fig.  2.  Lattice  parameters  of  Cu2Cd1_xZn*SnSe4  as  a  function  of  x.  The  formulas  show 
the  relationship  between  the  lattice  parameter  and  the  zinc  content  x,  where  the 
subscript  S  refers  to  the  stannite-type  structure  and  K  refers  to  the  kesterite-type 


the  kesterite-type  samples.  The  kesterite-type  samples  have 
shorter  Se-Sn  bonds  and  longer  Se-Cu  bonds  than  those  of  the 
stannite-type  samples. 

The  DSC  curves  for  Cu2CdxZnj_xSnSe4  presented  in  Fig.  4  reveals 
that  the  melting  point  (Tm)  increases  with  increasing  zinc  content  x 
for  both  the  stannite-  and  the  kesterite-type  compounds.  The  slight 
deviation  for  the  samples  with  x  =  0.55, 0.60,  and  0.65  may  be  due  to 
the  existence  of  atomic  disorder  in  the  structure  of  these  samples.  The 
melting  points  for  Cu2CdSnSe4  and  Cu2ZnSnSe4  are  1040  I<  and 
1061  K,  respectively,  which  agrees  with  the  values  reported  by 
Matsushita  et  al.  [20],  The  melting  point  of  Cu2-IMV-VI4  compounds 
has  been  determined  to  decrease  linearly  with  increasing  mean 
atomic  weight  (M)  as  follows:  Tm(°C)  =  -8.72 M  +  1500.  Our  results 
indicate  the  following  relationships:  Tms(°C)  =  — 4.2(4)M  +  850(5) 
and  TmK(0C)  —  -9.5(9  )M  + 1263(5),  where  the  subscripts  refers  to  the 
stannite-type  structure  and  I<  refers  to  the  kesterite-type  structure. 

3.2.  Thermoelectric  properties 

The  temperature  dependencies  of  the  electrical  transport  prop¬ 
erties  for  Cu2Cdi  xZnxSnSe4  are  shown  in  Fig.  5.  The  electrical 
conductivity  increases  for  the  samples  with  x  =  0  to  0.3  with 
increasing  temperature  but  decreases  slightly  for  the  others 


Fig.  3.  Atomic  distance  of  Cu2Cd,  „ZnJ(SnSe4  as  a  function  of  x. 
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(Fig.  5(a)).  The  electrical  conductivity  also  increases  with  the 
increasing  Zn  content  x,  except  for  fluctuations  in  the  cases  of  the 
samples  with  0.55  <  X  <  0.70  (Fig.  5(a),  inset).  The  electrical  con¬ 
ductivity  for  Cu2ZnSnSe4  in  this  work  was  1.53  x  104  S  m-1  at  300  K 
and  0.78  x  104  S  m-1  at  720  K,  which  are  much  higher  than  the 
values  reported  by  Liu  et  al.  [  1 2  ]  but  lower  than  those  reported  by  Shi 
et  al.  [  13  ] .  The  electrical  conductivity  for  the  Cu2CdSnSe4  in  this  work 
was  1.6  x  102  Sm  1  at  300 1<  and  9.1  x  102  S  m  1  at  720  K,  which  are 
slightly  lower  than  the  values  reported  by  Liu  et  al.  [12]  and  those  of 
the  samples  prepared  by  nano-powder  and  chemical  methods  [15]. 

The  differences  in  the  reported  electrical  transport  properties  may  U 
be  due  to  slight  differences  in  the  chemical  composition,  micro-  l> 
structure,  or  trace  second  phase  that  result  from  the  different  syn-  w 
thesis  methods.  The  electrical  conductivities  of  these  Cu-based 
quaternary  compounds  have  been  calculated  as  being  generated 
from  the  hybridization  of  Cu  3d  and  Se  4 p  orbitals  near  the  valence- 
band  maximum  (VBM)  [16],  The  Zn  atoms  doped  into  the 
Cu2CdSnSe4  structure  decrease  the  Cu-Se  atomic  distance,  which 
may  increase  the  concentration  of  carriers.  Flowever,  the  Cu2Se4 
tetrahedral  array  is  a  layer  in  the  plane  parallel  to  the  ab  plane  in  the 
structure  of  Cu2CdSnSe4  separated  by  insulating  CdSnSe2  tetrahe¬ 
dral  slabs.  In  contrast,  the  distribution  is  reticular  in  the  Cu2ZnSnSe4 
structure.  Disorder  of  Cu  and  Zn  atoms  in  the  kesterite-type  struc¬ 
ture  has  also  been  observed  [23,24],  which  may  be  due  to  the  close 
cationic  radii  of  Cu+  and  Zn2+.  When  the  Zn  atoms  are  doped  into  the 
structure  of  Cu2CdSnSe4,  the  Cu  and  Zn  ions  may  also  be  disordered, 
which  would  provide  more  conduction  channels,  thereby  resulting 
in  greater  electrical  conductivity.  Lattice  distortion  affects  the 
banding  length  and  thus  the  distribution  of  electrons  near  the  Fermi 
surface,  which  may  be  another  reason  for  the  relationship  between 
the  electrical  conductivity  and  the  zinc  content  x. 

The  Seebeck  coefficients  and  power  factors  increase  with 
increasing  temperature  for  all  of  the  samples  (Fig.  5(b  and  c)).  The 
positive  value  of  the  Seebeck  coefficient  indicates  P-type  conduc¬ 
tion  for  all  the  samples.  Associated  with  the  higher  electrical  con¬ 
ductivity  for  the  samples  with  x  =  0.55,  0.60,  or  0.65,  the  Seebeck 
coefficients  for  these  samples  are  relatively  lower  (Fig.  5(b),  inset). 

Except  for  these  samples,  the  Seebeck  coefficients  for  the  other 
samples  decrease  with  increasing  Zn  content  x.  The  Seebeck  co¬ 
efficients,  i.e.,  109  pV  K”1  at  300  K  and  219  pV  K-1  at  720  K  for 
Cu2ZnSnSe4  and  262  pV  K-1  at  300  K  and  349  pV  K"1  at  720  K  for 
Cu2CdSnSe4,  are  higher  than  those  of  some  well-known  thermo¬ 
electric  materials.  The  relationship  between  resistivity  and  the 


Seebeck  coefficient  agrees  with  the  traditional  relationship  be¬ 
tween  the  concentration  and  mobility  of  the  carrier.  The  power 
factor  (P)  also  increases  with  the  Zn  content  x,  except  for  fluctua¬ 
tions  in  the  cases  of  the  samples  with  0.55  <  X  <  0.70  (Fig.  5(c)). 


r(K) 


r(K) 


T(K) 


Fig.  5.  Temperature  dependence  of  the  electrical  conductivity  (a);  the  Seebeck  coef¬ 
ficient  (b);  and  the  power  factor  (c).  The  insets  show  the  electrical  conductivity,  the 
Seebeck  coefficient,  and  the  power  factor  as  functions  of  the  Zn  content  x  at  300  K  and 
720  K.  The  legends  for  (a)  and  (b)  are  the  same  as  those  for  (c). 


xr/W/mK)  k  (W/m/K) 


FS.  Liu  et  al.  /  Intermetallics  55  (2014)  15-21 


The  temperature  dependence  of  the  thermal  conductivity  of 
Cu2Cdi_xZnxSnSe4  is  demonstrated  in  Fig.  6(a  and  b).  The  com¬ 
pounds  with  the  stannite-type  structure  exhibit  lower  thermal 
conductivity  (Fig.  6(a))  compared  with  the  compounds  with  the 
kesterite-type  structure  (x  =  0.7,  0.8,  0.9, 1.0)  at  both  300  K  and 
720  K  (Fig.  6(b),  inset).  The  lattice  thermal  conductivity  (k;)  was 
calculated  by  subtracting  the  carrier  thermal  conductivity  ( kc )  from 
the  total  thermal  conductivity  (k).  The  carrier  thermal  conductivity 
kc  was  estimated  using  the  Wiedemann-Franz  relation  kc  =  LzaT, 
where  a  Lorenz  constant  of  Lz  —  2.0  x  10-8  W  Q  I<  2  was  applied  for 
the  estimation.  The  results  indicate  that  kc  is  very  small  for  all  of  the 
samples,  with  a  maximum  value  of  0.11  W  rrr1  K-1  in  Cu2ZnSnSe4 
at  720  K  (Fig.  6(c))  and  therefore  does  not  greatly  affect  the  total 
thermal  conductivity.  The  variation  of  the  lattice  thermal  conduc¬ 
tivity  with  the  Zn  content  x  exhibits  the  same  tendency  as  the  total 
thermal  conductivity  (Fig.  6(b),  inset).  This  result  indicates  that  the 
thermal  conductivity  is  mainly  dominated  by  the  lattice  conduc¬ 
tivity.  Therefore,  the  lower  thermal  conductivity  of  Cu2CdSnSe4  is 
attributed  to  its  greater  lattice  distortion  compared  to  that  in 
Cu2ZnSnSe4,  which  enhances  phonon  scattering  by  the  lattice. 

According  to  the  Keyes  expression  [25],  when  the  temperature  is 
greater  than  the  Debye  temperature,  the  lattice  thermal  conduc¬ 
tivity  (k/)  is  dominated  primarily  by  phonon-phonon  scattering: 

T  R?'2  Tm3/2p2/3  D  Tm3/2p%3 

K‘  3r2e3WV3  M7/ 6  MV  6  ’ 

where  Tm  is  the  melting  point,  M  is  the  mean  atomic  weight,  y  is  the 
Gruneisen  constant,  e  is  the  amplitude  of  interatomic  thermal  vi¬ 
bration,  R  is  the  ideal  gas  constant,  N0  is  Avogadro's  number,  and  p 
is  the  density.  Keyes  stated  that  B  can  be  regarded  as  a  constant  to 
estimate  the  lattice  thermal  conductivity  at  high  temperatures.  On 
the  basis  of  the  data  listed  in  Table  1,  the  calculated  value  of 
kq  =  Tm3VpVi /m7/6  increases  with  increasing  Zn  content  x,  which 
agrees  with  the  relationship  between  thermal  conductivity  and  the 
zinc  content  x  at  720  K  (Table  1  and  Fig.  6(b),  inset). 

The  ZT  calculated  from  the  data  above,  as  shown  in  Fig.  7,  in¬ 
creases  sharply  with  increasing  temperature  because  the  power 
factor  increases  while  the  thermal  conductivity  decreases.  The  ZT 
for  Cu2Cd1_xZnxSnSe4  increases  with  increasing  Zn  content  x  up  to 
x  =  0.60  and  then  fluctuates  with  a  further  increase  in  x.  Because  of 
its  relatively  higher  electrical  conductivity  and  lower  thermal 


T  (K) 


of  the  Zn 


F.S.  Liu  et  al.  /  Intermetallics  55  (2014)  15-21 


conductivity,  Cu2Cdo.4Zno.6SnSe4  has  a  maximum  ZT  value  of  0.23  at 
720  K  (Fig.  6,  inset).  The  ZT  values  of  the  solid  solutions  in  this 
system  primarily  depend  on  the  electrical  conductivity.  The  elec¬ 
trical  conductivity  increased  from  910  S  m-1  to  7812  S  m  ’, 
whereas  the  power  factor  increased  more  than  three-fold  as  the 
zinc  content  x  increased  from  0  to  1  at  720  K.  However,  the  thermal 
conductivity  only  increased  from  0.91  W  m-1  K-1  to 

1.32  W  m-1  K-1  at  the  same  temperature.  This  result  may  imply 
that  the  low  electrical  conductivity  leads  to  a  low  ZT  value  in  the 
solid  solution.  The  electrical  conductivity  of  Cu2.iZno.gSnSe4 
reached  33,200  S  m-1  by  Cu-doping,  and  the  ZT  value  reached  0.9  at 
860  I<  [12],  Therefore,  increasing  the  carrier  concentration  is  one  of 
the  most  effective  routes  to  further  enhance  the  thermoelectric 
performances  of  the  Cu2CdSnSe4  and  Cu2ZnSnSe4  solid  solutions. 

4.  Conclusions 

In  Cu2Cd|_xZnxSnSe4  solid  solutions,  the  lattice  parameters  for 
both  the  stannite-type  (x  =  0,  0.1,  0.2,  0.3,  0.5,  0.55,  0.6,  0.65)  and 
the  kesterite-type  (x  =  0.7,  0.8,  0.9, 1.0)  structures  linearly  evolve 
with  increasing  zinc  content  x.  Smaller  zinc  contents  result  in 
smaller  c/a  ratios  and  greater  lattice  distortion,  leading  to  lower 
thermal  conductivity.  The  Seebeck  coefficient  tends  to  decrease 
with  increasing  Zn  content,  whereas  the  electrical  and  thermal 
conductivities  tend  to  increase.  A  maximum  ZT  of  0.23  was  obtained 
for  the  sample  of  Cu2Cdo.4Zno.6SnSe4  at  720  K. 
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